Zebrafish as a model to investigate a biallelic gain-of-function variant in MSGN1, associated with a novel skeletal dysplasia syndrome.

BACKGROUND/OBJECTIVES: Rare genetic disorders causing specific congenital developmental abnormalities often manifest in single families. Investigation of disease-causing molecular features are most times lacking, although these investigations may open novel therapeutic options for patients. In this study, we aimed to identify the genetic cause in an Iranian patient with severe skeletal dysplasia and to model its molecular function in zebrafish embryos. RESULTS: The proband displays short stature and multiple skeletal abnormalities, including mesomelic dysplasia of the arms with complete humero-radio-ulna synostosis, arched clavicles, pelvic dysplasia, short and thin fibulae, proportionally short vertebrae, hyperlordosis and mild kyphosis. Exome sequencing of the patient revealed a novel homozygous c.374G > T, p.(Arg125Leu) missense variant in MSGN1 (NM_001105569). MSGN1, a basic-Helix-Loop-Helix transcription factor, plays a crucial role in formation of presomitic mesoderm progenitor cells/mesodermal stem cells during early developmental processes in vertebrates. Initial in vitro experiments show protein stability and correct intracellular localization of the novel variant in the nucleus and imply retained transcription factor function. To test the pathogenicity of the detected variant, we overexpressed wild-type and mutant msgn1 mRNA in zebrafish embryos and analyzed tbxta (T/brachyury/ntl). Overexpression of wild-type or mutant msgn1 mRNA significantly reduces tbxta expression in the tailbud compared to control embryos. Mutant msgn1 mRNA injected embryos depict a more severe effect, implying a gain-of-function mechanism. In vivo analysis on embryonic development was performed by clonal msgn1 overexpression in zebrafish embryos further demonstrated altered cell compartments in the presomitic mesoderm, notochord and pectoral fin buds. Detection of ectopic tbx6 and bmp2 expression in these embryos hint to affected downstream signals due to Msgn1 gain-of-function. CONCLUSION: In contrast to loss-of-function effects described in animal knockdown models, gain-of-function of MSGN1 explains the only mildly affected axial skeleton of the proband and rather normal vertebrae. In this context we observed notochord bending and potentially disruption of pectoral fin buds/upper extremity after overexpression of msgn1 in zebrafish embryos. The latter might result from Msgn1 function on mesenchymal stem cells or on chondrogenesis in these regions. In addition, we detected ectopic tbx6 and bmp2a expression after gain of Msgn1 function in zebrafish, which are interconnected to short stature, congenital scoliosis, limb shortening and prominent skeletal malformations in patients. Our findings highlight a rare, so far undescribed skeletal dysplasia syndrome associated with a gain-of-function mutation in MSGN1 and hint to its molecular downstream effectors.

Spatiotemporal control of pattern formation during somitogenesis.

Spatiotemporal patterns widely occur in biological, chemical, and physical systems. Particularly, embryonic development displays a diverse gamut of repetitive patterns established in many tissues and organs. Branching treelike structures in lungs, kidneys, livers, pancreases, and mammary glands as well as digits and bones in appendages, teeth, and palates are just a few examples. A fascinating instance of repetitive patterning is the sequential segmentation of the primary body axis, which is conserved in all vertebrates and many arthropods and annelids. In these species, the body axis elongates at the posterior end of the embryo containing an unsegmented tissue. Meanwhile, segments sequentially bud off from the anterior end of the unsegmented tissue, laying down an exquisite repetitive pattern and creating a segmented body plan. In vertebrates, the paraxial mesoderm is sequentially divided into somites. In this review, we will discuss the most prominent models, the most puzzling experimental data, and outstanding questions in vertebrate somite segmentation.

Stochastic gene expression and environmental stressors trigger variable somite segmentation phenotypes.

Mutations of several genes cause incomplete penetrance and variable expressivity of phenotypes, which are usually attributed to modifier genes or gene-environment interactions. Here, we show stochastic gene expression underlies the variability of somite segmentation defects in embryos mutant for segmentation clock genes her1 or her7. Phenotypic strength is further augmented by low temperature and hypoxia. By performing live imaging of the segmentation clock reporters, we further show that groups of cells with higher oscillation amplitudes successfully form somites while those with lower amplitudes fail to do so. In unfavorable environments, the number of cycles with high amplitude oscillations and the number of successful segmentations proportionally decrease. These results suggest that individual oscillation cycles stochastically fail to pass a threshold amplitude, resulting in segmentation defects in mutants. Our quantitative methodology is adaptable to investigate variable phenotypes of mutant genes in different tissues.

The Role of Fibroblast Growth Factor Signaling in Somitogenesis.

Fibroblast growth factor (FGF) signaling is conserved from cnidaria to mammals (Ornitz and Itoh, 2022) and it regulates several critical processes such as differentiation, proliferation, apoptosis, cell migration, and embryonic development. One pivotal process FGF signaling controls is the division of vertebrate paraxial mesoderm into repeated segmented units called somites (i.e., somitogenesis). Somite segmentation occurs periodically and sequentially in a head-to-tail manner, and lays down the plan for compartmentalized development of the vertebrate body axis (Gomez et al., 2008). These somites later give rise to vertebrae, tendons, and skeletal muscle. Somite segments form sequentially from the anterior end of the presomitic mesoderm (PSM). The periodicity of somite segmentation is conferred by the segmentation clock, comprising oscillatory expression of Hairy and enhancer-of-split (Her/Hes) genes in the PSM. The positional information for somite boundaries is instructed by the double phosphorylated extracellular signal-regulated kinase (ppERK) gradient, which is the relevant readout of FGF signaling during somitogenesis (Sawada et al., 2001; Delfini et al., 2005; Simsek and Ozbudak, 2018; Simsek et al., 2023). In this review, we summarize the crosstalk between the segmentation clock and FGF/ppERK gradient and discuss how that leads to periodic somite boundary formation. We also draw attention to outstanding questions regarding the interconnected roles of the segmentation clock and ppERK gradient, and close with suggested future directions of study.

A design logic for sequential segmentation across organisms.

Multitudes of organisms display metameric compartmentalization of their body plan. Segmentation of these compartments happens sequentially in diverse phyla. In several sequentially segmenting species, periodically active molecular clocks and signaling gradients have been found. The clocks are proposed to control the timing of segmentation, while the gradients are proposed to instruct the positions of segment boundaries. However, the identity of the clock and gradient molecules differs across species. Furthermore, sequential segmentation of a basal chordate, Amphioxus, continues at late stages when the small tail bud cell population cannot establish long-range signaling gradients. Thus, it remains to be explained how a conserved morphological trait (i.e., sequential segmentation) is achieved by using different molecules or molecules with different spatial profiles. Here, we first focus on sequential segmentation of somites in vertebrate embryos and then draw parallels with other species. Thereafter, we propose a candidate design principle that has the potential to answer this puzzling question.

Human stem cell models unravel mechanisms of somite segmentation.

In vitro models to study human somitogenesis, the formation of the segmented body plan, have so far been limited.1 Two papers in Nature now report the creation of pluripotent stem cell (PSC)-derived 3D culture systems that recapitulate the formation of somite-like structures and help gain insights into this developmental process.2,3.

Using single-molecule fluorescence in situ hybridization and immunohistochemistry to count RNA molecules in single cells in zebrafish embryos.

Taming gene expression variability is critical for robust pattern formation during embryonic development. Here, we describe an optimized protocol for single-molecule fluorescence in situ hybridization and immunohistochemistry in zebrafish embryos. We detail how to count segmentation clock RNAs and calculate their variability among neighboring cells. This approach is easily adaptable to count RNA numbers of any gene and calculate transcriptional variability among neighboring cells in diverse biological settings. For complete details on the use and execution of this protocol, please refer to Keskin et al. (2018),1 Zinani et al. (2021),2 and Zinani et al. (2022).3.

Periodic inhibition of Erk activity drives sequential somite segmentation.

Sequential segmentation creates modular body plans of diverse metazoan embryos1-4. Somitogenesis establishes the segmental pattern of the vertebrate body axis. A molecular segmentation clock in the presomitic mesoderm sets the pace of somite formation4. However, how cells are primed to form a segment boundary at a specific location remains unclear. Here we developed precise reporters for the clock and double-phosphorylated Erk (ppErk) gradient in zebrafish. We show that the Her1-Her7 oscillator drives segmental commitment by periodically lowering ppErk, therefore projecting its oscillation onto the ppErk gradient. Pulsatile inhibition of the ppErk gradient can fully substitute for the role of the clock, and kinematic clock waves are dispensable for sequential segmentation. The clock functions upstream of ppErk, which in turn enables neighbouring cells to discretely establish somite boundaries in zebrafish5. Molecularly divergent clocks and morphogen gradients were identified in sequentially segmenting species3,4,6-8. Our findings imply that versatile clocks may establish sequential segmentation in diverse species provided that they inhibit gradients.

Patterning principles of morphogen gradients.

Metazoan embryos develop from a single cell into three-dimensional structured organisms while groups of genetically identical cells attain specialized identities. Cells of the developing embryo both create and accurately interpret morphogen gradients to determine their positions and make specific decisions in response. Here, we first cover intellectual roots of morphogen and positional information concepts. Focusing on animal embryos, we then provide a review of current understanding on how morphogen gradients are established and how their spans are controlled. Lastly, we cover how gradients evolve in time and space during development, and how they encode information to control patterning. In sum, we provide a list of patterning principles for morphogen gradients and review recent advances in quantitative methodologies elucidating information provided by morphogens.

Gene copy number and negative feedback differentially regulate transcriptional variability of segmentation clock genes.

Timely progression of a genetic program is critical for embryonic development. However, gene expression involves inevitable fluctuations in biochemical reactions leading to substantial cell-to-cell variability (gene expression noise). One of the important questions in developmental biology is how pattern formation is reproducibly executed despite these unavoidable fluctuations in gene expression. Here, we studied the transcriptional variability of two paired zebrafish segmentation clock genes (her1 and her7) in multiple genetic backgrounds. Segmentation clock genes establish an oscillating self-regulatory system, presenting a challenging yet beautiful system in studying control of transcription variability. In this study, we found that a negative feedback loop established by the Her1 and Her7 proteins minimizes uncorrelated variability whereas gene copy number affects variability of both RNAs in a similar manner (correlated variability). We anticipate that these findings will help analyze the precision of other natural clocks and inspire the ideas for engineering precise synthetic clocks in tissue engineering.

Regulatory mechanisms ensuring coordinated expression of functionally related genes.

Coordinated spatiotemporal expression of large sets of genes is required for the development and homeostasis of organisms. To achieve this goal, organisms use myriad strategies where they form operons, utilize bidirectional promoters, cluster genes, share enhancers among genes by DNA looping, and form topologically associated domains and transcriptional condensates. Coexpression achieved by these different strategies is hypothesized to have functional importance in minimizing gene expression variability, establishing dosage balance to ensure stoichiometry of protein complexes, and minimizing accumulation of toxic intermediate metabolites. By combining gene-editing tools with computational modeling, recent studies tested the advantages of adjacent genes located in pairs and clusters. We propose that with the advancement of gene editing, single-cell sequencing, and imaging tools, one could readily test the functional importance of different coexpression strategies in a variety of biological processes.

A 3-D Tail Explant Culture to Study Vertebrate Segmentation in Zebrafish.

Vertebrate embryos pattern their major body axis as repetitive somites, the precursors of vertebrae, muscle, and skin. Somites progressively segment from the presomitic mesoderm (PSM) as the tail end of the embryo elongates posteriorly. Somites form with regular periodicity and scale in size. Zebrafish is a popular model organism as it is genetically tractable and has transparent embryos that allow for live imaging. Nevertheless, during somitogenesis, fish embryos are wrapped around a large, rounding yolk. This geometry limits live imaging of PSM tissue in zebrafish embryos, particularly at higher resolutions that require a close objective working distance. Here, we present a flattened 3-D tissue culture method for live imaging of zebrafish tail explants. Tail explants mimic intact embryos by displaying a proportional slowdown of axis elongation and shortening of rostrocaudal somite lengths. We are further able to stall axis elongation speed through explant culture. This, for the first time, enables us to untangle the chemical input of signaling gradients from the mechanistic input of axial elongation. In future studies, this method can be combined with a microfluidic setup to allow time-controlled pharmaceutical perturbations or screening of vertebrate segmentation without any drug penetration concerns.

Pairing of segmentation clock genes drives robust pattern formation.

Gene expression is an inherently stochastic process1,2; however, organismal development and homeostasis require cells to coordinate the spatiotemporal expression of large sets of genes. In metazoans, pairs of co-expressed genes often reside in the same chromosomal neighbourhood, with gene pairs representing 10 to 50% of all genes, depending on the species3-6. Because shared upstream regulators can ensure correlated gene expression, the selective advantage of maintaining adjacent gene pairs remains unknown6. Here, using two linked zebrafish segmentation clock genes, her1 and her7, and combining single-cell transcript counting, genetic engineering, real-time imaging and computational modelling, we show that gene pairing boosts correlated transcription and provides phenotypic robustness for the formation of developmental patterns. Our results demonstrate that the prevention of gene pairing disrupts oscillations and segmentation, and the linkage of her1 and her7 is essential for the development of the body axis in zebrafish embryos. We predict that gene pairing may be similarly advantageous in other organisms, and our findings could lead to the engineering of precise synthetic clocks in embryos and organoids.

Regulatory Network of the Scoliosis-Associated Genes Establishes Rostrocaudal Patterning of Somites in Zebrafish.

Gene regulatory networks govern pattern formation and differentiation during embryonic development. Segmentation of somites, precursors of the vertebral column among other tissues, is jointly controlled by temporal signals from the segmentation clock and spatial signals from morphogen gradients. To explore how these temporal and spatial signals are integrated, we combined time-controlled genetic perturbation experiments with computational modeling to reconstruct the core segmentation network in zebrafish. We found that Mesp family transcription factors link the temporal information of the segmentation clock with the spatial action of the fibroblast growth factor signaling gradient to establish rostrocaudal (head to tail) polarity of segmented somites. We further showed that cells gradually commit to patterning by the action of different genes at different spatiotemporal positions. Our study provides a blueprint of the zebrafish segmentation network, which includes evolutionarily conserved genes that are associated with the birth defect congenital scoliosis in humans.

Spatial Fold Change of FGF Signaling Encodes Positional Information for Segmental Determination in Zebrafish.

Signal gradients encode instructive information for numerous decision-making processes during embryonic development. A striking example of precise, scalable tissue-level patterning is the segmentation of somites-the precursors of the vertebral column-during which the fibroblast growth factor (FGF), Wnt, and retinoic acid (RA) pathways establish spatial gradients. Despite decades of studies proposing roles for all three pathways, the dynamic feature of these gradients that encodes instructive information determining segment sizes remained elusive. We developed a non-elongating tail explant system, integrated quantitative measurements with computational modeling, and tested alternative models to show that positional information is encoded solely by spatial fold change (SFC) in FGF signal output. Neighboring cells measure SFC to accurately position the determination front and thus determine segment size. The SFC model successfully recapitulates results of spatiotemporal perturbation experiments on both explants and intact embryos, and it shows that Wnt signaling acts permissively upstream of FGF signaling and that RA gradient is dispensable.

Noise in the Vertebrate Segmentation Clock Is Boosted by Time Delays but Tamed by Notch Signaling.

Taming cell-to-cell variability in gene expression is critical for precise pattern formation during embryonic development. To investigate the source and buffering mechanism of expression variability, we studied a biological clock, the vertebrate segmentation clock, controlling the precise spatiotemporal patterning of the vertebral column. By counting single transcripts of segmentation clock genes in zebrafish, we show that clock genes have low RNA amplitudes and expression variability is primarily driven by gene extrinsic sources, which is suppressed by Notch signaling. We further show that expression noise surprisingly increases from the posterior progenitor zone to the anterior segmentation and differentiation zone. Our computational model reproduces the spatial noise profile by incorporating spatially increasing time delays in gene expression. Our results, suggesting that expression variability is controlled by the balance of time delays and cell signaling in a vertebrate tissue, will shed light on the accuracy of natural clocks in multi-cellular systems and inspire engineering of robust synthetic oscillators.

Reduced dosage of ß-catenin provides significant rescue of cardiac outflow tract anomalies in a Tbx1 conditional null mouse model of 22q11.2 deletion syndrome.

The 22q11.2 deletion syndrome (22q11.2DS; velo-cardio-facial syndrome; DiGeorge syndrome) is a congenital anomaly disorder in which haploinsufficiency of TBX1, encoding a T-box transcription factor, is the major candidate for cardiac outflow tract (OFT) malformations. Inactivation of Tbx1 in the anterior heart field (AHF) mesoderm in the mouse results in premature expression of pro-differentiation genes and a persistent truncus arteriosus (PTA) in which septation does not form between the aorta and pulmonary trunk. Canonical Wnt/ß-catenin has major roles in cardiac OFT development that may act upstream of Tbx1. Consistent with an antagonistic relationship, we found the opposite gene expression changes occurred in the AHF in ß-catenin loss of function embryos compared to Tbx1 loss of function embryos, providing an opportunity to test for genetic rescue. When both alleles of Tbx1 and one allele of ß-catenin were inactivated in the Mef2c-AHF-Cre domain, 61% of them (n = 34) showed partial or complete rescue of the PTA defect. Upregulated genes that were oppositely changed in expression in individual mutant embryos were normalized in significantly rescued embryos. Further, ß-catenin was increased in expression when Tbx1 was inactivated, suggesting that there may be a negative feedback loop between canonical Wnt and Tbx1 in the AHF to allow the formation of the OFT. We suggest that alteration of this balance may contribute to variable expressivity in 22q11.2DS.

Spatial gradients of protein-level time delays set the pace of the traveling segmentation clock waves.

The vertebrate segmentation clock is a gene expression oscillator controlling rhythmic segmentation of the vertebral column during embryonic development. The period of oscillations becomes longer as cells are displaced along the posterior to anterior axis, which results in traveling waves of clock gene expression sweeping in the unsegmented tissue. Although various hypotheses necessitating the inclusion of additional regulatory genes into the core clock network at different spatial locations have been proposed, the mechanism underlying traveling waves has remained elusive. Here, we combined molecular-level computational modeling and quantitative experimentation to solve this puzzle. Our model predicts the existence of an increasing gradient of gene expression time delays along the posterior to anterior direction to recapitulate spatiotemporal profiles of the traveling segmentation clock waves in different genetic backgrounds in zebrafish. We validated this prediction by measuring an increased time delay of oscillatory Her1 protein production along the unsegmented tissue. Our results refuted the need for spatial expansion of the core feedback loop to explain the occurrence of traveling waves. Spatial regulation of gene expression time delays is a novel way of creating dynamic patterns; this is the first report demonstrating such a control mechanism in any tissue and future investigations will explore the presence of analogous examples in other biological systems.

Short-lived Her proteins drive robust synchronized oscillations in the zebrafish segmentation clock.

Oscillations are prevalent in natural systems. A gene expression oscillator, called the segmentation clock, controls segmentation of precursors of the vertebral column. Genes belonging to the Hes/her family encode the only conserved oscillating genes in all analyzed vertebrate species. Hes/Her proteins form dimers and negatively autoregulate their own transcription. Here, we developed a stochastic two-dimensional multicellular computational model to elucidate how the dynamics, i.e. period, amplitude and synchronization, of the segmentation clock are regulated. We performed parameter searches to demonstrate that autoregulatory negative-feedback loops of the redundant repressor Her dimers can generate synchronized gene expression oscillations in wild-type embryos and reproduce the dynamics of the segmentation oscillator in different mutant conditions. Our model also predicts that synchronized oscillations can be robustly generated as long as the half-lives of the repressor dimers are shorter than 6 minutes. We validated this prediction by measuring, for the first time, the half-life of Her7 protein as 3.5 minutes. These results demonstrate the importance of building biologically realistic stochastic models to test biological models more stringently and make predictions for future experimental studies.

Cited3 activates Mef2c to control muscle cell differentiation and survival.

Vertebrate muscle development occurs through sequential differentiation of cells residing in somitic mesoderm - a process that is largely governed by transcriptional regulators. Our recent spatiotemporal microarray study in zebrafish has identified functionally uncharacterized transcriptional regulators that are expressed at the initial stages of myogenesis. cited3 is one such novel gene encoding a transcriptional coactivator, which is expressed in the precursors of oxidative slow-twitch myofibers. Our experiments placed cited3 into a gene regulatory network, where it acts downstream of Hedgehog signaling and myoD/myf5 but upstream of mef2c. Knockdown of expression of cited3 by antisense morpholino oligonucleotides impaired muscle cell differentiation and growth, caused muscle cell death and eventually led to total immotility. Transplantation experiments demonstrated that Cited3 cell-autonomously activates the expression of mef2c in slow myofibers, while it non-cell-autonomously regulates expression of structural genes in fast myofibers. Restoring expression of cited3 or mef2c rescued all the cited3 loss-of-function phenotypes. Protein truncation experiments revealed the functional necessity of C-terminally conserved domain of Cited3, which is known to mediate interactions of Cited-family proteins with histone acetylases. Our findings demonstrate that Cited3 is a critical transcriptional coactivator functioning during muscle differentiation and its absence leads to defects in terminal differentiation and survival of muscle cells.